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Abstract. We analyze the introduction of dissipative effects in the study of the 
dynamical Casimir effect. We consider a toy model for an electromagnetic cavity that 
contains a semiconducting thin shell, which is irradiated with short laser pulses in 
£h ' order to produce periodic oscillations of its conductivity. The coupling between the 

quantum field in the cavity and the microscopic degrees of freedom of the shell induces 
q-i dissipation and noise in the dynamics of the field. We argue that the photon creation 

process should be described in terms of a damped oscillator with nonlocal dissipation 
and colored noise. 
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q ■ 1. Introduction 
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The motion of a neutral body distorts the quantum electromagnetic vacuum and may 
induce photon creation. This fact was first pointed out many years ago p. The radiation 
of a single mirror is, in general, an extremely small effect. For this reason, during 
the seventies, the radiation of moving mirrors was mainly analyzed because its formal 



relation with black hole evaporation [2j. However, it was later pointed out that, in a 
resonant cavity, photon creation could be enhanced by secular effects: if the length of 
the cavity oscillates with a frequency which is twice the frequency of an eigenmode of 
the static cavity, resonant effects may produce a number of photons that, under ideal 
conditions, grows exponentially with time [3]. 

Although intensively studied from the theoretical point of view, this "motion 
induced radiation" is extremely difficult to measure. On the one hand, cavities with 
a very high Q are needed in order to have an exponential growth during a significative 
amount of time. On the other hand, and more importantly, for microwave cavities 
the mirror should oscillate at extremely high frequencies, of the order of Ghz. This is 
a challenge from an experimental point of view, although there are feasible proposals 
based on the use of nanoresonators in order to accomplish it [I]. 

Photon creation can occur whenever there are time dependent external conditions: 
moving boundaries and/or time dependent electromagnetic properties. These broad 
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class of phenomena is usually denominated "Dynamical Casimir Effect" (DCE) (for a 
complete list of publications see j5]). In particular, a setup that has attracted both 
theoretical and experimental attention is the possibility of using short laser pulses in 
order to produce periodic variations of the conductivity of a semiconductor layer placed 
inside a microwave cavity. The fast changes in the conductivity induce a periodic 
variation in the effective length of the cavity, and therefore the creation of photon 
pairs [6]. This setup has been analyzed at the theoretical level [8], and there is an 
ongoing experiment aimed at the detection of the motion induced radiation [9]. More 
recently, the possibility of changing the effective length of a superconducting coplanar 
waveguide terminated by a SQUID using a time dependent magnetic field has also been 
put forward [TO] . 

In this paper we will be mainly concerned with the irradiated-semiconductor setup. 
It has been pointed out that dissipative effects may play an important role in this 
experiment, that could induce a significant deviation from the expected exponential 
growth in the number of created photons [8] . To our knowledge, up to now there is no 
theoretical model that includes, from first principles, effects of dissipation and noise. 
There are, on the other hand, phenomenological models pj] to estimate these effects. 
The idea is the following: in a resonant situation, if the spectrum of the cavity is not 
equidistant, the photons are created in a single mode. The dynamics of this mode 
can be described, in the absence of dissipation, in terms of a harmonic oscillator with 
time dependent frequency. A large number of photons is created due to parametric 
resonance. A phenomenological way of introducing dissipation and noise is to describe 
the temporal evolution of the modes in terms of a damped oscillator. Mathematical 
consistency for quantization requires the introduction of an additional stochastic force, 
the noise. Physically, this is nothing but the well known fluctuation-dissipation relation. 

The aim of this paper is to provide a first step towards a description of the 
DCE from first principles, including dissipation and noise. The natural arena for this 
description is the theory of quantum open systems. We will consider the electromagnetic 
field as our "system", while the degrees of freedom on the semiconductor will be the 
"environment" . Under very general assumptions, is it possible to show that the dynamics 
of the electromagnetic field will be described by a Langevin equation (Section 3). In 
general, this equation will include a nonlocal dissipative term and a colored noise, both 
concentrated on the position of the layer. We will use this Langevin equation to describe 
the dynamics of the resonant mode (Section 4). We will see that the evolution of a single 
mode is in general more complex than that of a damped oscillator with white noise, and 
that all the effects come from the non trivial boundary conditions at the position of the 
slab. 

In order to avoid unnecessary complications we will consider several simplifying 
assumptions. On the one hand we will consider a quantum scalar field in 1 + 1 dimensions 
instead of the electromagnetic field in a 3 + 1 cavity. On the other hand, we will assume 
that this quantum field is linearly coupled to the degrees of freedom of the layer. Finally, 
we will consider a very thin semiconductor layer. Even with these simplifications, we 
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hope the resulting model to contain the main physics of the problem. This model can 
be considered as a generalization of the theory described in Ref. |7], reviewed in Section 
2, to the dissipative case. 

2. A model without dissipation 

In this section we will review the model of Ref. [7], which we will use as starting point 
to include dissipative effects in the rest of the paper. 

We consider a massless scalar field within a cavity with perfect conducting walls. 
For simplicity, we consider the 1+1 dimensional case, where we consider a cavity of size 
L. At the mid point of the cavity (x = L/2) a thin film of semiconducting material is 
located. We model the conductivity properties of such material by a potential V{t): The 
ideal limit of perfect conductivity corresponds to V — > oo, and V — > to a 'transparent' 
material. This potential varies between a minimum value, Vo, and a maximum Vm ax . 
The Lagrangian of the scalar field within the cavity is given by 

£ = \d^d^ - ^-8(x - L/2)<f\ (1) 

where 5(x) is the one-dimensional Dirac delta function. The use of an infinitely thin 
film is justified as long as the width of the slab is much smaller than the wavelengths of 
the relevant modes in the cavity. The corresponding Lagrange equation reads, 

{dl - d 2 t )<P = V(t)S(x - L/2)4>. (2) 

As shown in Ref. [12]. this model describes plane-polarized electromagnetic field 
propagating normally to an infinitesimally thin jellium-type plasma sheet. The strength 
of the potential is given by 

V = Ann s e 2 /m , (3) 

where n s is the surface charge density of free carriers in the sheet, e is the electron 
charge, and m is the effective mass of the free carriers. When a laser field suddenly 
impinges on the plasma sheet, it produces time- dependent changes in the surface charge 
density n s (t) that induce a time variation in the conductivity V(t). 

For the sake of clarity we divide the cavity into two regions: region I (0 < x < L/2) 
and region II (L/2 < x < L). Perfect conductivity at the edges of the cavity imposes 
the following boundary conditions 

= 0,t) = (j)u(x = L,t) = 0. (4) 

The presence of the semiconducting film introduces a discontinuity in the spatial 
derivative of the field, while the field itself remains continuous, 

<k(x = L/2,t) = Mx = L/2,t), 

d x Ux = L/2, t) - d x <h.{x = L/2, t) = -V{t)<f>{x = L/2, t). (5) 

This can be seen by integrating out the field equation (Eq. (j2J)) in the neighborhood of 
the film. 
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We consider the set of functions 

i(j m (x,t) = sin (k m (t)x) , (6) 
where k m (t) is the m-th positive solution to the following transcendental equation 

2fc m = -y(t)tan(%^). (7) 



2 

Note that ip m depends on t through k m (t). 
Let us define 

m[ , ) \ -<iP m (x-L,t) L/2<x<L W 

These functions satisfy the boundary conditions Eq.(jl]) and Eq.([5]), and the 
orthogonality relations (9 m ,9 n ) = [1 — sm{k m {t)L) /kmij^L] 5 m , n , where we have used 
the usual inner product in the interval [0,L]. 

For t < the slab is not irradiated, consequently V is independent of time and 
has the value Vq. The modes of the quantum scalar field that satisfy the Klein Gordon 
equation (T2J) are 

g— iu> m t 

u m (x,t) = -==V m (x,0) , (9) 

V ^ m 

where u m = k m , and k m is the m-th solution to Eq.([7]) for V = V . At t = the 
potential starts to change in time and the set {k m } of the eigenfrequencies of the cavity 
acquires a time dependence through Eq.([7]). 
We expand the field operator (p as 

cf)(x, t) = Y^ [b m u m (x, t) + b m u m (x, t)] , (10) 
in 

where b m and b m are annihilation and creation operators, respectively. There is another 
set of solutions of the Klein Gordon equation that satisfy all boundary conditions and 
have a node at x = L/2. For this reason their dynamics is not affected by the presence 
of the slab. Moreover, these modes are decoupled from the modes u m thanks to the 
orthogonality conditions. Therefore, we will only consider the evolution of the modes 
given in Eq.Q. 

For t > we write the expansion of the field mode u s as 

u s (x,t>o) = ^p r { : ) m m (^t)- (ii) 

m 

Replacing this expression into — d^)u s = we find 

P^ + kn{t)Pn^ = ~ ^ [(^P^&jn + ^m^m) 9 ml + Pm^mdmn > (12) 

m 

where the coefficients gmn read 



9mn ~ (*„, *„) V dkl ' * n ' • 
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Note that u m = k m (0). 

We are interested in the number of photons created inside the cavity. Hence we 
focus in resonance effects induced by periodic oscillations in the conductivity V(t), which 
translates into effective periodic changes in the modes of the scalar field. Therefore we 
start by considering a time dependent conductivity given by 

V(t) = V + (V max - V Q ) f(t) , (14) 

where f{t) is a periodic and non negative function, f{t) = f(t + T) > 0, that vanishes 
at t = and attains its maximum at f{r e ) = 1. In each period, f(t) describes the 
excitation and relaxation of the semiconductor produced by the laser pulse. Typically, 
the characteristic time of excitation r e is the smallest time scale and satisfies r e <C T. 
We expand f(t) in a Fourier series 

oo 

f(t) = f + J2fj™<jM + c j ), (15) 

J=l 

where Q = 2ir/T. Since r e is the smallest time scale, on general grounds we expect the 
first T/r e terms in the above series to be relevant. 

Under certain constraints, large changes in V induce only small variations in k 
through the transcendental relation between k and V (see Eq.([7j)). In this case, a 
perturbative treatment is valid and a linearization of such relation is appropriate. 
Accordingly we write 

k n (t) = k° n (l + e n f(t)), (16) 

where e n is obtained after replacing Eqs. (I14II16P into Eq. ([7]) and expanding it to first 
order in e n . The result is 

The restriction for the validity of the perturbative treatment is V L ^> V max /V > 1. 
These conditions are satisfied for realistic values of L, Vq, and Vm ax . Indeed, from 
Eq.fl3]), and using known values for the conductivities of good conductors, we can fix 
Vmax = 10 16 m _1 . When the laser field is not applied we can set Vq — 10 10 m _1 — 10 13 m _1 , 
the range of values for different semiconductors. For a cavity of size L x ~ 10~ 2 m, and 
when the ratio between the maximum and minimum conductivities is in the range 
10 3 < V max /V < 10 6 , we obtain from Eq.flE]) small values of e, 10" 8 < e < 10" 2 . It is 
worth noticing that we are interested in low eigenfrequencies, for which k(t) ~ 0(L~ l ). 
Nonetheless the perturbative treatment is also valid for k ~ 0(V). 

In what follows we will only consider expressions to first order in e n . To analyze 
the possibility of parametric resonance we write the time dependent frequency given in 
Eq. (DED as, 

k n (t) = ~k° n (l + e n (f-fo)), (18) 
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where = + e n /o) is a 'renormalized' frequency. The equation for the coefficients 
Pm\t) (Eq. (fl2]) ) can now we written to first order in e n as, 

P^ + (k° n M s) = - 2e n (k° n f(f - / )pM 

[2P^e m k°J + PW^/] ^ + 0(e 2 ). (19) 



This equation describes a set of coupled harmonic oscillators with periodic frequencies 
and couplings. It is of the same form as the equations that describe the modes of a scalar 
field in a three dimensional cavity with an oscillating boundary. A naive perturbative 
solution of previous equations in powers of e n breaks down after a short amount of time 
(this happens for particular values of the external frequency such that there is a resonant 
coupling with eigenfrequencies of the cavity). In order to analyze the long time behavior 
of the solutions the equations can be solved using multiple scale analysis (MSA) [T3] . 

As shown in detail in Ref.[7], if the spectrum of the cavity is not equidistant, then 
the photons are created mainly in the resonant mode. This means that this particular 
mode decouples from the rest, and the subsequent quantum dynamics can be described 
starting from the equation 

Pi n) + (k° n ) 2 P^ + 2e n (k° n ) 2 (f - fo)Pi n) = , (20) 

that is, a harmonic oscillator with time dependent frequency. When Q = 2k®, the 
number of created photons grows exponentially due to parametric resonance. 



3. Quantum open systems: The Langevin equation for the field 

In the previous section, the starting point was the action given in Eq.([T]), which describes 
a quantum field with a time-dependent mass term concentrated on the position of the 
slab. In a more realistic model, the massless scalar field should be coupled to the 
microscopic degrees of freedom on the slab, and its dynamics will be described in terms 
of the effective action that results after integrating the microscopic degrees of freedom. 

In static situations, it is enough to compute the in-out or the Euclidean effective 
action, and to read the zero point energy of the system from the generating functional. 
However, if we are interested in the temporal evolution of the quantum field, it is 
necessary to compute the in-in or CTP effective action [TJ]. In this direction, we will 
consider the scalar field as a quantum open system, which is coupled to an external 
environment composed by the internal degrees of freedom in the slab. Therefore, we 
will consider model composed by a classical action given by 

S[<f>, q n ] = S [<f>] + S [q n ] + S int [(f), q n ], (21) 

where So[</>] and S [q n ] are the free actions for the massless scalar field and for the infinite 
set of harmonic oscillators with which we model the microscopic degrees of freedom on 
the slab, 

So[<t>] = \j ^xd^dy, (22) 
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n 

and the interaction term is assumed as 

S- m t[<P,q n ] = ~ / d 2 x\(t)<j)(t,x)q n (t)5(x - -), (24) 

n 

where \(t) is the time-dependent coupling between the scalar field and the microscopic 
degrees of freedom on the slab. With this time dependent coupling we are modelling the 
time dependent conductivity of the slab, which could be produced by a time dependent 
charge carrier density within the semiconducting slab. 

Our system of interest is the massless scalar field. Therefore, we will proceed 
integrating out the harmonic oscillator variables q n in order to obtain an effective 
description of the evolution of the field. This non-unitary evolution is described, mainly, 
by two environmentally induced effects: dissipation and noise. After tracing over the 
environment oscillators, we are able to write an effective action for the quantum open 
system [15] of the form 

A[<j>, 4>>] = S [cf>} - SoW] + SA[4>, 0'], (25) 

where and are each of the CTP branches for the field. The influence (Feynmann- 
Vernon) action is given by 



f )£(*',*') 

^)A(t',x') (26) 

where we have defined A(t,x) = </>(t,x) — <f)'(t,x) and E(i, x) = l/2[<f)(t,x) + (f>'(t,x)]. 

The influence action is in general non-local and complex, and contains all the 
information about the environment (i.e. temperature, spectral density, dissipation 
rates, and noise properties). The kernel in the real part of 5 A is the dissipation kernel 
D(t;t'), and N(t;t') in the imaginary part is related with the quantum fluctuations or 
noise coming from the external bath. These kernels are usually related by means of a 
fluctuation-dissipation relation. The time- dependent coupling strenght is included in 
the definition of the noise and dissipation kernels in Eq. fj26j) . On general grounds, one 
can expect that D(t; f ) = X(t)X(t')D(t - t') and N(t; t') = X(t)X(t')N(t - f ). 

Different environments can be considered. Each of them is characterized by a 
spectral density. The ohmic environment is the most studied case in the literature and 
produces a dissipative force that is proportional to the velocity. The supraohmic case, on 
the one hand, is generally used to model the interaction between defects and phonons in 
metals and also to mimic the interaction between a charge and its own electromagnetic 
field [IE]. On the other hand, the quantum behaviour of free electrons in mesoscopic 
systems is affected by their interaction with the environment, which, for example in such 
cases, consist of other electrons, phonons, photons or scatterers. Which environment 



5A[<P,<P'}= J d 2 x J d 2 x'A(t,x)D(t;t')5(x- ^)5{x' 



L. 



+ i d 2 x d 2 x'A(t,x)N(t;t')5(x - -)S(x' 
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is more relevant for the dissipative phenomena generally depends on the temperature. 
For instance, the temperature dependence of the weak-localization correction to the 
conductivity reveals in metals that electron-electron interactions dominate over the 
phonon contribution to decoherence at the low temperature regime. Therefore, in order 
to mimic a real material for the slab in the cavity, it will be relevant to consider spectral 
densities such the ones mentioned above. Only in the particular case in which the 
coupling is constant X(t) = A , and we consider an ohmic enviroment at very high 
temperature, the noise kernel is N(t — t') ~ Aq5(£ — t') and the dissipation kernel reads 
D(t — t') ~ Aq5(£ — t'). Therefore, the influence action is local in time. It is worth to 
remark that, for these models, the fluctuation-dissipation relation reads 

/+oo 
dsD(s), (27) 
-oo 

which is the classical Einstein formula. 

In the general case, the equation of motion of the scalar field can be obtained from 
the effective action A[(p, by |^|^ = ^/ = 0: 

□0 + 5(x ~^)J dsD(t; s) = 0. (28) 

This equation corresponds to an average over all the possible noise realizations. In order 
to extract the noise information from the effective action, it is possible to write down 
the imaginary part of the influence action in terms of a stochastic noise force J 7 , which 
is coupled to the main system of interest, and it is defined by means off a Gaussian 
probability distribution PfJ 7 ] like 

P[JT] = Njrexpl~J ds J ds'F(s)N{s; s')- 1 ^') | , (29) 

where Njr is a normalization constant. Therefore, the imaginary part of 5 A can be 
written in terms of the stochastic force source as 



= exp^jy s ^ds'A(s,^)N(s;s')A(s',^)y (30) 

Finally, from Eq.( l30]) we can derive an equation of motion for the massless scalar 
field that takes into account the noise. This is a Langevin-like equation, where the full 
dynamics of the field is also determined by the presence of the stochasticity induced by 
the environment via J 7 . The Langevin equation now reads as 

□0 + 8{x - \) ^ dsD{t- S )0(|, s) = JF{t)8{x - |). (31) 

The stochastic noise T is characterised by the probability distribution -PfJ 7 ], and 
also by the correlation functions 



<^(t)> = ; {T{f)T{f)) = N(t;t>) 
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4. The dynamics of the resonant mode 

The aim of this section is to obtain a dynamical equation for the evolution of a single 
mode of the quantum field, starting from the effective equation derived in the previous 
section. We have already seen that, in the absence of dissipation, if the conductivity of 
the slab has periodic changes with frequency Q, only the resonant mode with k = ^fl 
is relevant at long times. We will now assume that this is true even in the presence 
of dissipation. Although reasonable, this is an unjustified assumption, and we hope to 
further analyze it in a forthcoming publication. Even with this simplification, we will 
see that the derivation of the dynamical equation of the mode is a non-trivial task. 
Let us write the dissipation kernel as 

D(t; s) = V S(t -s) + d{t; s) . (32) 

In the particular case d(t; s) = (V(t) — Vo)5(t — s), with V(t) a periodic function of 
frequency Q we recover the model of Section 2. In that case, we have shown that even 
when V(t) differs several orders of magnitude from Vq, the oscillations in the wavenumber 
have a relative amplitude e„ C 1 (see the discussion after Eq.(|T71)). Therefore, we will 
assume that the main contribution to the wavenumber comes from the term proportional 
to Vq, and we will treat the contributions of both d(t; s) and the noise kernel as small 
corrections. It is worth to emphasize that we are not assuming that D(t, s) ~ VoS(t — s), 
but that the effect of d(t; s) and the noise on the instantaneous wavenumber k(t) can 
be treated perturbatively. 

The resonant mode of the field can be written as 

u(x,t) = P(t)V(x,t) (33) 

where \l/(x, t) is given in Eq.(|6]). In order to simplify the notation, from now on we omit 
the subindex m, since we are considering a single mode. Inserting this particular mode 
into the Langevin equation for the field, we find that the boundary condition at the 
position of the slab becomes 

2P(t)k(t) cos[^] = ~f da D(t, s)P(s) sin[^] + yJ±F{t) . (34) 

It is worth to stress that, when the dissipation kernel is local and in the absence of noise, 
the amplitude of the mode P(t) factorizes and the boundary condition fixes the value 
of the time dependent wavenumber k(t). In the general case, k(t) becomes a nonlocal 
function of P{t). To see this explicitly, we solve Eq.( l34"l) assuming that d(t; s) and ^(t) 
induce small time dependent corrections to k(t). We write k(t) = k° + Ak(t), where k° 
is the solution of the unperturbed (static) problem, i.e. 

2k° = -VbtanS . (35) 
To first order in d(t; s) and T we have 



k 



o 



Ak(t)P(t) = 

Vo[1 + m + ii^k] 



[ dsd(t;s)P(s) - F(t) 
Jo 



(36) 
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where J- = \[L]2 cosec(/c°L/2). As anticipated, the nonlocal boundary condition 
becomes a nonlocal relation between Ak(t) and P(t). 

The next step is to insert the field mode into the Klein Gordon equation of the 
field, at both sides of the slab. Of course the single mode will not be an exact solution 
of the Klein Gordon equation, which would involve the coupling to the other modes. 
Therefore, in order to get the dynamical equation for P(t), we impose (*&(x, t), Ou) = 0. 
After some simple calculations we get, 

P + ((A; ) 2 + 2k°Ak{t))P + g {A) {2PAk + PAk) = (37) 

where corresponds to gmm is defined in Section 2, and can be evaluated to lowest 
order, setting d(t; s) — and T = 0. Note that in this approximation this coupling 
constant is independent of time. 

Let us analyze the resulting equation for the field mode. To begin with, it is 
important to note that the term proportional to Ak(t)P(t) contains nonlocal dissipation 
an additive noise. Moreover, if we write d(t; s) = (V(t) — V )5(t — s) + d(t\ s) the 
local term corresponds to a time dependent frequency. Depending on the details of 
the environment, we expect that, at high temperatures, d(t — s) ~ 5(t — s) so the 
usual dissipative term proportional to P is recovered. Moreover, in this limit the 
fluctuation-dissipation theorem implies a white noise. In any other case, we expect 
nonlocal dissipation and colored noise. 

In the absence of dissipation, the last term in the lhs of Eq. (I37j) . proportional to 
g( A \ can be neglected when the external frequency satisfies f2 = 2k°. However, in 
general this is not the case. Noting that 2PAk + PAk = (PAk) — PAk, it is simple 
to find an expression for this additional term by taking two derivatives of Eq.()36l). The 
result is a rather complicated expression: 

k° 



2PAk + PAk = -PAk 



V [l + ( fc0 ) 2L i 



V 



+ d{t;t)P+ [ dsd$d(t;s)P(s)-F(t) 
Jo 



P(d(t;t)+dtd(t;t)) 



(38 



The main consequence is that, after inserting this result into Eq. (1371) . the final equation 
for P(t) can be written as 

P + oo 2 cS (t)P + f ds d eS (t, s)P(s) = F cS (t) . (39) 
Jo 

In summary, we have shown that the amplitude of the mode satisfies a Langevin 
equation with nonlocal dissipation and colored noise. The dissipation (d e g) and noise 
(J- e s) kernels are related to the ones in the Langevin equation for the field derived in 
Section 3, although they are not exactly the same because of the presence of the term 
proportional to g^ . All the information about dissipation and noise comes from the 
nonlocal boundary condition that the field satisfies on the slab. 
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5. Conclusions 

In this paper we have argued that the natural approach to analyze the dynamical Casimir 
effect is that of the quantum open systems. The coupling between the degrees of freedom 
of the vacuum field and those of the (imperfect) mirrors generate an infuence functional 
for the vacuum field that contains the information about dissipation and noise. The 
dynamics of the field is described by a Langevin equation with dissipative and noise 
kernels concentrated on the position of the mirror. 

When the system is under the influence of periodic, time dependent external 
conditions, it is possible to have parametric amplification in some of the modes of 
the field. If the spectrum of the unperturbed system is not equidistant, one expects 
the amplification to occur for a single mode. Assuming that this is the case even in 
the presence of dissipation, we have shown that the dynamics of the mode can also be 
described by a Langevin equation with nonlocal dissipation and colored noise. This 
equation is similar to the one derived in the context of quantum brownian motion, but 
now the frequency and the (nonlocal) kernels have a periodic time dependence. 

We have considered a number of simplifications in order to illustrate the above 
points: we worked with a scalar field in 1 + 1 dimensions, and we assumed that 
the parametric amplification can be described in terms of a single mode. The last 
assumption deserve further investigation. Moreover, we have not attempted to compute 
the influence of dissipation and noise on the number of photons created. This seems 
to be a rather difficult task, because of the complexity of the Langevin equation for 
the mode. Some additional simplifications could be necessary in order to estimate these 
effects (along the lines of Ref.JTTJ, for instance). However we think that we have clarified, 
from a conceptual point of view, the origin of the description of the DCE in terms of a 
noisy and damped oscillator. 

This work has been supported by Universidad de Buenos Aires, CONICET and 
ANPCyT. 
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